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Abstract 

We analyze the origin of the quasiclassical realm from the no-boundary proposal for the universe's 
quantum state in a class of minisuperspace models. The models assume homogeneous, isotropic, 
closed spacetime geometries, a single scalar field moving in a quadratic potential, and a fundamental 
cosmological constant. The allowed classical histories and their probabilities are calculated to 
leading semiclassical order. We find that for the most realistic range of parameters analyzed a 
minimum amount of scalar field is required, if there is any at all, in order for the universe to 
behave classically at late times. If the classical late time histories are extended back, they may 
be singular or bounce at a finite radius. The ensemble of classical histories is time symmetric 
although individual histories are generally not. The no-boundary proposal selects inflationary 
histories, but the measure on the classical solutions it provides is heavily biased towards small 
amounts of inflation. However, the probability for a large number of efoldings is enhanced by the 
volume factor needed to obtain the probability for what we observe in our past light cone, given 
our present age. Our results emphasize that it is the quantum state of the universe that determines 
whether or not it exhibits a quasiclassical realm and what histories are possible or probable within 
that realm. 
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I. INTRODUCTION 



The inference is inescapable from the physics of the last eighty years that we live in a 
quantum mechanical universe. If so, the universe has a quantum state. A theory of that state 
is as important a challenge for fundamental physics as a theory of the dynamics. Providing 
that theory, and testing its observational predictions, are the goals of quantum cosmology. 

A central prediction of the universe's quantum state is the classical spacetime that is a 
manifest fact of the present universe. Predicting classical spacetime is a constraint on the 
theory of the state because we can no more expect to find classical predictions following from 
a general state in quantum gravity than we can in the non-relativistic quantum mechanics 
of a particle. Histories exhibit classical correlations in time only when they are suitably 
coarse-grained and then only for particular kinds of states (e.g 

The probabilities for the alternative classical histories of a quantum universe answer 
questions such as the following: Is approximate homogeneity and isotropy likely? Is the 
probability high for sufficient inflation to explain the present spatial flatness? Is a homo- 
geneous thermodynamic arrow of time likely? What is the probability that the universe 
bounced at a minimum radius above the Planck scale in the past? 

This paper is concerned with the classical histories predicted by the no-boundary wave 
function of the universe (NBWF) [2| in homogeneous, isotropic minisuperspace models with 
a fundamental cosmological constant and a single scalar field moving in a quadratic potential. 
Many of our results and conclusions together with speculations concerning their extensions 
to other models have been summarized in This paper presents the detailed derivations 
of these and deals exclusively with quadratic potentials for the scalar field. 

By way of introduction we now briefiy sketch the standard procedure for classical predic- 
tion in quantum cosmology. More details will be found in Section II, and derivations in the 
context of generalized quantum theory in 0]. 

States in quantum cosmology are represented by wave functions on the superspace of 
three-geometries and spatial matter field configurations. For the homogeneous, isotropic, 
spatially closed, minisuperspace models with one scalar field that are the subject of this pa- 
per, wave functions depend on the scale factor b determining the size of the spatial geometry 
and the value x of the homogeneous scalar field. Thus, \E' = x). 

A class of states of particular interest are those whose wave functions can be approximated 
to leading order in h in some region of superspace by the semiclassical form (or superpositions 
of such forms) 

^{b, x) ~ exp{[-/«(6, x) + ^S{b, x)]/h} (1.1) 

with both and 5* real. When 5* varies rapidly and Ir varies slowly such wave functions 
predict an ensemble of suitably coarse-grained Lorentzian histories with high probabilities 
for correlations in time governed by classical deterministic laws for spacetime geometry and 
the matter field. This requirement on the gradients of Ir and S is called the classicality 
condition. When it is satisfied the action 5* determines the ensemble as in familiar Hamilton- 
Jacobi theory. Classical histories not contained in the ensemble have zero probability in this 
approximation. The classical histories that are members of the ensemble have probabili- 
ties proportional to exp[— 27/^(6, x)]//i]. In this way particular states can predict classical 
spacetime. 

The no-boundary wave function for these models is defined by the sum-over-histories 

^{b,x) = [ 6a6(Pexp{-I[a{T),^{r)]/h). (1.2) 
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Here, a(r) and 0(r) are the histories of the scale factor and matter field and /[a(r), 0(t)] is 
their Euclidean action. The sum is over cosmological geometries that are regular on a disk 
with only one boundary at which a(r) and 0(r) take the values b and x- The integration is 
carried out along a suitable complex contour C which ensures the convergence of fll.2p and 
the reality of the result [5|. 

For some ranges of b and x it niay happen that the integral in fll.2p can be approximated 
by the method of steepest descents. Then the wave function will be well approximated by 
a sum of terms of the form (11.11) — one for each extremizing history (a(r), 4>{t)) matching 
(6, x) on the boundary of the manifold and regular elsewhere. In general these solutions 
will be complex — "fuzzy instantons". For each contribution x) is the real part of 
the action /[a(r),0(r)] evaluated at the extremizing history and —S{b,x) is the imaginary 
part. When the classicality condition is satisfied the no-boundary wave function predicts an 
ensemble of classical histories as described above. 

Two key points should be noted about this prescription for classical prediction: (1) The 
NBWF provides probabilities for entire classical histories. It therefore supplies a classical 
history measure, i.e. a measure on classical phase space that is conserved along the classical 
trajectories. (2) The histories in the classical ensemble are not the same as the extremizing 
histories that provide the steepest descents approximation to the integral (11. 2p defining the 
NBWF. 

We apply this prescription for classical prediction to the NBWF in minisuperspace models 
with a cosmological constant A and a quadratic potential of the form V{^) = (l/2)m^$^. 
Our main aim is determining the ensemble of classical cosmologies predicted by the no- 
boundary proposal and the probabilities of its members. From these we calculate the 
probabilities for whether universe bounces or is singular, for whether it expands forever 
or recollapses, for the magnitude of any time asymmetries, for different amounts of matter 
content, and for different amounts of inflation. 

Our detailed conclusions are given in Section [IX] , but were it necessary to single out 
just two they would be the following: (1) Not all classical behaviors of the universe are 
allowed by the no-boundary proposal and for some ranges of model parameters no classical 
behavior is predicted at all. The manifest existence of the quasiclassical realm in this 
universe is therefore an important, non-trivial, constraint on theories of its initial quantum 
state. (2) The probability for significant inflation depends sensitively on the limitations of 
the classical ensemble arising from the classicality condition and also on the limited scale of 
our observations in a large universe and the values of cosmological parameters such as the 
present age. 

The predictions of the no-boundary quantum state were extensively analyzed in minisu- 
perspace models in the '80s and '90s (see e.g. We are unable to give anything like a 
complete survey of this work, but relevant papers include the following: Classicality condi- 
tions were studied in similar models in 0,0, B especially in connection with the amount 
of inflation predicted by the NWBF. For related work on this question see jo], T3] and [ll[ 
for the influence of higher order quantum corrections which are neglected here. Complex 
solutions were studied in detail by [l2|, 13, 14]; this paper relies heavily on these works. 

What is new is the following: (1) A better understanding of the prescription for classical 
prediction attained by providing a firmer foundation through the probabilities for histories 
provided by generalized quantum theory 0, [l5|. (2) A more complete analysis of the com- 
plex solutions that provide the semiclassical approximation to the no-boundary proposal 
and their connection to the probabilities for classical cosmologies. (3) The inclusion of a 
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fundamental cosmological constant which both generahzes the discussion and simphfies it. 
(4) A derivation of the probabihties for what we observe within our past hght cone from the 
probabihties for the ensemble of entire classical histories predicted by the NBWF. 

This paper is organized as follows: In Section [ITl we review the prescription for extracting 
the predictions for classical cosmologies from a wave function of the universe. Section IIIII 
describes the no-boundary wave function and its semiclassical approximation. The homo- 
genenous, isotropic minisuperspace models that are the focus of this paper are laid out in 
detail in Section I1V[ Section |V] analyzes the complex fuzzy instantons that provide the 
semiclassical approximation to the NBWF. The predicted ensemble of classical ensemble 
and the probabilities of its members are discussed in Section |Vll Section IVIII discusses con- 
ditional (top down) probabilities relevant for the amount of inflation of the universe given 
our present observations. Section IVIIll discusses the arrow of time in those histories that 
bounce at a minimum radius in our past. Section [IX] contains our conclusions. 



II. CLASSICAL PREDICTION IN QUANTUM COSMOLOGY 

A quantum system behaves classically when the probability is high for histories of its mo- 
tion that exhibit patterns of classical correlation in time governed by deterministic dynamical 
laws. The Moon can be said to move on a classical orbit when the quantum mechanical 
probability is high for histories of coarse grained positions of the Moon's center of mass that 
obey Newton's laws. The universe behaves classically when the quantum probability is high 
for histories of coarse grained geometry and matter fields that are correlated in time by the 
Einstein equation. 

In this section we summarize the prescription for classical prediction in quantum cosmol- 
ogy. The context is the minisuperspace models with homogeneous, isotropic geometries and 
homogeneneous scalar field studied in this paper. These rules can be derived in a quantum 
framework that predicts probabilities for sets of alternative, coarse-grained histories of cos- 



mological geometries and matter fields whether or not they behave classically [15|, [IGj, \11 
Alternatively the rules can be motivated as a simple extension of the analogous algorithm 
derived in non-relativistic quantum mechanics. For brevity here we defer both of these 
arguments to a separate paper [3]. 

A detailed discussion of our minisuperspace models will be given in Section IIV[ But for 
the present discussion it is sufficient to note that they lie in the class specified by a classical 
action S of the form: 



S[N{X),q^{X)] = K / dXN 



1 /I dq^\ f 1 dq 



2 \NdXj\Nd\' 



(2.1) 



Here, the q^ are a set of coordinates for the minisuperspace of homogeneous, isotropic 
three-geometries and homogeneous three-dimensional field configurations. For our models 
the scale factor b and the field value x ci^re the q^. iV is a multiplier ensuring reprametriza- 
tion invariance. (The hat is for consistency with later notation.) Histories are curves in 
minisuperspace specified by giving these coordinates as a function of a parameter A, viz 
q^{\)- The metric^ on superspace GAsiQ^) and the potential V{q^) specify the model. The 
constant K is fixed by scaling conventions for the variables and the conventions for S. 



Up to factors this is the inverse DeWitt metric. 
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The action (12. ip is invariant under reparametrizations of the histories q^{\)- As a con- 
sequence there is a constraint relating the coordinates and their conjugate momenta pa- 
This can be found by varying (12. ip with respect to A^(A) and expressing the result in terms 
of and pa- The result can be put in the form: 

H{pA, g^) = ^G^^'PAPB + V(g^) = 0. (2.2) 

In quantum cosmology the state of the universe is specified by giving a wave function on 
superspace. For minisuperspace models this is "^{q^). In this paper that is the no-boundary 
wave function (II. 2p . All wave functions satisfy an operator implementation of the classical 
constraint (12. 2p 

H (-^^^, ^(/) = (-f + V(g^)) ^(q^) = 0. (2.3) 

(We retain the factors of h for later convenience in discussing classicality.) This is the 
Wheeler-DeWitt equation for these models. There is a conserved current associated with 
the Wheeler-DeWitt equation 

This will play an important role in defining the probabilities for histories. 

Suppose that in some region of superspace the wave function of the universe has the 
approximate semiclassical form (or is a sum of such forms) 

^(/) ~ A(g^)e±'^(''^)/'* (2.5) 

where S{q^)/h varies rapidly over the region and A{q^) varies slowly. Under these cir- 
cumstances the Wheeler-DeWitt equation (12. 3p requires that S{q'^) satisfies the classical 
Hamilton- Jacobi equation to a good approximation 

In a suitable coarse-graining the only histories that have significant probability are the 
classical histories corresponding to the integral curves of S{q'^)/h (e.g. [4]). These are curves 
q^{X) which satisfy 

In short, wave functions of the form (12. 5p with rapidly varying S{q°')/h and slowly varying 
A(g"^) predict an ensemble of classical Lorentzian cosmological histories. 

Consider any surface in minisuperspace that is spacelike with respect to the metric Gab 
and has unit normal ua- We assume that the relative probability density p of classical 
histories passing through this surface is the component of the conserved current (12. 4p along 
the normal if it is positive. In leading order in h this is 

p{q^)^J.n=\A{q^)\'Vr.S{q^) (2.8) 

in any region in which it is positive. The first order in h implications of the Wheeler-DeWitt 
equation (12. 3p ensure that these probabilities are constant along classical trajectories. Thus 



6 



the formula (12. 8p could be evaluated on any spacelike surface with the same result for the 
probabilities of the classical histories that intersect it.^ 

Several key points should be noted about this prescription for classical prediction: 

• The no-boundary wave function provides probabilities for entire four-dimensional clas- 
sical histories. When the strong energy condition is not satisfied (as for the present 
models) these may bounce at a minimum radius. If that radius is large enough we 
expect a classical extrapolation from present data to be a good approximation over 
the whole history of the universe — from an infinite volume in the past to an infi- 
nite volume in the future. Alternatively when extrapolated from present data some 
classical histories assigned probabilities may have initial or final singularities or both. 

• Singularities in the extrapolation of classical solutions do not signal the break down 
of quantum mechanics but rather of the approximation. In particular the NBWF 
predicts probabilities for the classical description of late time observables such as 
CMB fluctuations whatever happens to an extrapolation of that classical description. 
That is because the NBWF predicts probabilities for histories not their initial data. 
In this sense the NBWF resolves classical singularities. 

• The histories in the classical ensemble are not the same as the extremizing histories that 
provide the steepest descents approximation to the integral (II. 2p defining the NBWF. 
The classical histories are real and Lorentzian. The extrema are generally complex — 
neither Euclidean nor Lorentzian except in very special cases. The classical histories 
may contract from an infinity in the past and reexpand to another one in the future; 
the no-boundary extrema can have only one infinity. Indeed, in the no-boundary case 
the classical histories and extremizing histories are on different manifolds. This clean 
separation into real classical histories and complex extremizing ones helps to clarify 
the meaning of both, and resolves issues that arise from their identification such as 
those dicsussed in [8|. 

From this point of view, a wave function of the universe is best thought of, not as an initial 
condition, but rather in a four dimensional sense as giving probabilistic weight to the possible 
four-dimensional histories of a quantum universe. 



III. THE NO-BOUNDARY WAVE FUNCTION AND ITS SEMICLASSICAL AP- 
PROXIMATIONS 

This section considers the steepest descents approximation to the NBWF for homogeneous 
isotropic minisuperspace models. It describes when this leads to a semiclassical form like 
(12. 5p which predicts an ensemble of classical Lorentzian histories with probabilities for each. 



^ This expression for probability of classical histories has been advocated by other authors, see e.g. |18l.ll' 
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A. Steepest descents approximation 



The NBWF is defined by a path integral over homogeneous field configurations and 
homogeneous isotropic metrics of the form 

ds^ = (3/A) [N\X)dX^ + a^{X)dnl] . (3.1) 

Here, dfl^ is the round metric on the unit three-sphere and the factor in front is a convenient 
normalization. The defining path integral has the specific form [cf. (11.21) ] 



■^{b,x) = '^iqn = J 5N6xexp{-I[N{X),x''{X)]/h) (3.2) 

where x^{X) = (a(A),0(A)) are histories of the scale factor and scalar field. The integral is 
over all (a;'^(A), A^(A)) that define regular geometries on a disk which match the values of 
= (6, x) on its boundary. The functional / is the Euclidean action which for our class of 
models is [cf. flO) ] 



I[N{X),x''{X)] =K / dXN 







(3.3) 



where the parameter values labeling the endpoints have been conventionally chosen to be 
and 1. The integration measure in (13.21) contains the usual apparatus of gauge fixing 
terms and their associated determinants made necessary by reparametrization invariance. 
We have left all of this unspecified because it will not be important in the leading steepest 
descents approximation. From now on in this paper \E'(g^) should be understood to be this 
no-boundary wave function. 

The steepest descents approximation to the path integral (13. 2p defining the NBWF 
starts with those paths that extremize the action in the class integrated over. Such paths 
(A^ext(A), x4t('^)) are solutions the equations of motion 

61 61 , , 

0, — -— = (3.4) 



6x^{X) ' 6N{X) 

that are regular on the disk and match the values on its boundary The explicit form of 
the equations (13. 4p will be displayed in Section HVl 

The steepest descents approximation to the NBWF is then given by 

vi/(g^)^5^exp[-Axt(/)A] (3.5) 

ext 

where the sum is over all extrema that contribute to the integral. The exponent has an 
expansion in powers of h of the form 

Axt(g^) = /ext(g^) + + ■ • • • (3.6) 

The leading order in this expansion is the Euclidean action evaluated at the extremum: 

/ext(/)=/[iVext(A),xi,(A)]. (3.7) 
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Next order corrections in h include terms like 

- {l/2)Tr\og{6^I) (3.8) 

where 6^1 is the operator resulting from the second variation of the action^. Factors arising 
from the measure would also contribute at this order. 

The Hamiltonian-Jacobi equation for the Euclidean action is the order consequence 
of the Wheeler-DeWitt equation (12.31) . The order implication is the conservation of 
the probabilities (12. 8p . We therefore should retain both orders in the steepest descents 
approximation to be consistent with these features. Traditonally the order h contributions 
are written as a prefactor to the exponential. Thus we have to write for the contribution to 
the wave function of one extremum 

*cxt(/) ~ Pext(/) exp[-/ext(g^)//^]. (3.9) 

From now on we will consider the extrema one at a time and drop the subscript "ext" 
that distinguished one from the other. 



B. Classicality 

There is no reason to assume that the leading steepest descents approximation will be 
given by a real path (iV(A), x"^(A)). The reality of the NBWF only means that the extrema 
must come in complex conjugate pairs. (See Section llll DI below for more on the consequences 
of this.) Indeed, we will show in Section IIII CI how the extremizing paths are necessarily 
complex. The action at an extremum will therefore have both real and imaginary parts 
which we write: 

/(/)=//?(/) -^5(g^)- (3.10) 

The first two orders in steepest descents approximation to the wave function therefore take 
the form 

^{q^) = A(g^)e'^(''^)/^ (3.11) 

with A given by 

A{q^) = P(g^)e-^«('?^)/^ (3.12) 

As reviewed in Section |Tll a wave function with the semiclassical form (13. lip in some 
region of minisuperspace predicts an ensemble of classical trajectories provided that S{q^)/h 
is rapidly varying and A{q^) is slowly varying. Assuming that P{q^) is slowly varying, a 
necessary condition for this is that the gradient of lR{q^) be small compared to the gradient 
of S{q^) in the coordinates that enter into the defining path integral and for which we 
expect to have classical equations of motion e.g {b,x) [4]. That is, 

\VaIr\ < |Va^|. (3.13) 
This is the classicality condition which plays a central role in our work. 



^ This correction leads to a well known prefactor which can be written as the inverse square root of the 
determinant of S^I 
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Further, as we saw in Section [TTl when the action S^q"^) satisfies the classical Hamilton- 
Jacobi equation 

^(VS)2 + V(g^)=0, (3.14a) 

the Lorentzian histories in the classical ensemble are the integral curves of S^q"^). Specif- 
ically, choosing N = i so the metrics (13.11) are Lorentizian, the integral curves obey the 
equations of motion 

PA{q^) = Gab^ = VA5(g^). (3.14b) 

Defined in (I3.10p as (minus) the imaginary part of the Euclidean action evaluated at an 
extremum, there is no reason to believe that eqs (I3.14p hold in all regions of minisuperspace, 
and indeed they do not. Rather the Euclidean action satisfies its own 'Euclidean Hamilton- 
Jacobi' equation 

-i(V/)^ + V(g^) = 0, (3.15a) 



with its own equation of motion 



1 flT^ 

Gabj^^ = VAl{q\ (3.15b) 



These follow directly from the definition (13.31) . 

Using (13.101) the real and imaginary parts of the Euclidean Hamilton- Jacobi relation 
(I3.15al) can be written as 

-\{VlRf + \{VSf + V(g^) = 0, (3.16a) 

V/r-V5 = 0. (3.16b) 

Eq (I3.16al) shows that the Hamilton- Jacobi equation (I3.14al) for S{q^) follows from the 
Euclidean Hamiltonian- Jacobi equation when the classicality condition (I3.13p holds because 
this implies 

\{yiR?\ « |(V5)2|. (3.17) 

But (I3.17P is not enough to guarantee classicality. It is sufficient for the Lorentizian 
Hamilton- Jacobi equation (I3.14ap . But it is not necessarily enough to ensure that the Eu- 
clidean equations of motion (]3.15bp reduce to their Lorentzian forms (]3.14bp . For that the 
stronger condition (I3.13P is necessary. 

The other consequence of the Euclidean Hamilton- Jacobi equation is (]3.16bp . This implies 



that iR^q^) is constant along the integral curves of S{q^). That is, each classical Lorentzian 
history is associated with a value of Ir. Indeed, in a two-dimensional minisuperspace like 
that of this paper, fl3.16bp implies that curves of constant Ir are integral curves of 5*. 

In principle the probability for any history can be calculated from the wave function 
^(g"^) without a semiclassical approximation. The semiclassical form is only a sufficient 
criterion for classicality. We will assume that once classical histories have been identified 
in a region of minisuperspace where the classicality condition holds they may be extended 
to regions where it does not hold using the classical equations of motion unless they become 
classically singular. It is plausible, for instance, that a bouncing universe whose radius never 
falls below the Planck length will remain classical throughout its history even if it can only 
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be identified by a steepest descents approximation in some regions of minisuperspace. That 
is an assumption wliicli can in principle be cliecked in tlie full quantum mechanical theory. 

We next turn to the probabilities predicted by the NBWF for the individual histories in 
the classical ensemble. According to (12.81) . the relative probability density p for classical 
histories passing through a spacelike surface in minisuperspace with unit normal is given 
by 

p(g^) = |P(g^)|V2^«(''")V„^ (3.18) 

where this expression is positive. For this to be a probability for histories it must be constant 
along the integral curves of S. As discussed in Section [Til the order approximation to the 
Wheeler-DeWitt equation (12. 3p ensures this. But since Ir is already constant along classical 
trajectories the rest of the measure (I3.18P must be separately constant. 

In this paper we will consider only the lowest order semiclassical approximation to the 
probabilities and ignore the prefactor P which arises in the next order. That is we discuss 
only the exp(— 2/^) contribution to the probabilities which is possible because it is conserved 
along classical histories. This 'approximation' is forced on us by our limited ability to 
compute the corrections to the leading term at this time. Higher order corrections can be 
crucial for some questions (e.g [ll|) but for the diagnosis of classicality we expect that the 
this lowest approximation is enough. 

C. Extrema of the Euclidean Action and Fuzzy Instantons 

We now return to a more detailed examination of the conditions determining the complex 
paths that extremize the action and the differential equations (13.41) which are necessary con- 
ditions for an extremum. The set of equations (13.41) consists of two second order differential 
equations for a(A) and 0(A) together with a constraint involving only their first deriva- 
tives. There are thus four real second order differential equations and two real constraint 
conditions. 

We first show that there are no free parameters in the boundary conditions that determine 
the solutions to these equations. The domain on which the equations are to be solved ranges 
from the center of symmetry of the geometry on the 3-disk at A = (called the 'South Pole' 
(SP)) to the boundary where b and x are specified at A = 1. The conditions for the geometry 
and field to be regular at the SP are 

a(0) = 0, 0'(O) = (3.19a) 

where a prime denotes a derivative with respect to A. The conditions at the boundary are 

a(l)=6, 0(1) =x (3.19b) 

where b and x a-re real. Eqs (I3.19P constitute eight real conditions for the four real second 
order equations for a(A) and 0(A). 

With a suitable choice of parametrization the multiplier A^(A) can be taken to be a 
complex constant N. For each N solve the second order differential equations with the 
boundary conditions (I3.19P for a(A) and 0(A). Then find the real and imaginary parts of 

so that the two real constraint conditions are satisfied. Were the equations linear the 
solutions would be determined. There may be more than one solution to this non-linear 
set of equations with the boundary conditions (I3.19P but there are no free parameters to 
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specify. Hence whether the solutions are real or complex is not up to us; it is determined 
by the equations and the boundary conditions. 

In the case of a cosmological constant and no scalar field there is a real solution consisting 
of a real Euclidean instanton with the geometry of half a round 4-sphere joined smoothly 
onto de Sitter space through a sphere of minimum radius [2()]. But as shown by a number of 
authors (see e.g. and, as will be verified here, there are no real extrema when the scalar 
field is non-zero except for special "false-vacuum" potentials. In general the extrema are 
necessarily complex. There is thus generally no meaningful notion of a Euclidean instanton 
nucleating the universe. But, as we will see in Section |Vl for a wide range of parameters 
the imaginary parts of extremizing solutions that satisfy the classicality condition fl3.13p 
are small. In this regime we can therefore think of these extremizing geometries as "fuzzy 
instantons" in which there is a transition from a real Euclidean geometry at the South Pole 
to an asymptotically real Lorentzian geometry at large volume. The transition is not sharp 
as in the zero scalar field case, but rather spread out over a region. We will give explicit 
examples in Section |Vl 

We stress again, however, that the complex fuzzy instantons that provide the semiclassi- 
cal approximation to the NBWF are distinct from the Lorentzian histories in the classical 
ensemble for which they provide probabilities through the real part of their complex action. 

D. Time Symmetry of the Classical Ensemble 

The metric Gab{(1^) and the potential V(g'^) that define the Euclidean action (13.31) are 
real analytic functions of their arguments for the models we will consider. So, therefore, are 
the coefficients in the equations (13. 4p which are the necessary conditions for an extremum 
of the action. The boundary conditions for their solutions (I3.19P are real. 

Therefore, for every extremum (A^(A), a(A), 0(A)), there is also a complex conjugate ex- 
tremum (Ar*(A),a*(A),0*(A)). If 

I{b,x) = IR{b,x)-^S{b,x) (3.20a) 
is the action for the first solution then the action for the second will be 

I{b,x)=IR{b,x)+^S{b,x)■ (3.20b) 

The real part of both actions is the same. Both extrema therefore count equally in their 
contribution to the steepest descents approximation to the NBWF. This shows explicitly 
that the NBWF is real in the semiclassical approximation. 

But the opposite signs for S in (I3.20p means that the momenta of a classical history 
passing through q"^ will be opposite in the two cases [cf. (12.70 ]. The classical ensembles of 
each extremum consist of histories that are time reversals of one another. Both will have the 
same probability because 1^(6, x) is the same for both. The individual histories need not 
be time symmetric [l3| and indeed we will find that they are not [cf . Figure [H] . But the 
ensemble of predicted classical histories is time symmetric in the sense that for any history 
in it, its time reversed is also a member with the same probability. 

E. Measures on Classical Phase Space 

A classical history measure is any function on phase space that is conserved along classical 
trajectories because of the classical equations of motion. Classical history measures on the 
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phase space of classical minisuperspace models have been adroitly employed by Gibbons and 
Turok [10| to analyze the probability of inflation in the absence of a theory of the universe's 
state. 

The predictions of the NBWF for an ensemble of classical histories provides a history 
measure on the classical phase space of minisuperspace models. The condition (13.14b|) be- 
tween the coordinates q"^ and the momenta pa shows the the NBWF measure is concentrated 
on a surface in classical phase space of half its dimension. One could think of the NBWF 
measure as a 5-function on this slice through phase space that assigns zero probability to 
those points not on it. 

As we will see in Section |Vl the NBWF surface does not slice through the whole of clas- 
sical phase space. Quantum mechanics assigns probabilities generally to decoherent sets of 
alternative histories of the universe. But only in special circumstances are the probabilities 
high for the correlations in time that define classical histories. A classical history therefore 
cannot be expected to pass through every point q^. The classicality condition (13.131) gen- 
erally specifies a boundary to the surface in phase space on which points corresponding to 
classical histories lie. 

This restriction of the ensemble of possible classical histories to a surface with boundary 
in phase space is already a powerful prediction of the NBWF whatever relative probabilities 
are predicted for histories in it. Eq. (I3.18P gives the predictions of the NBWF for the 
probabilities of the histories within the surface in phase space defining the classical ensemble. 
These probabilities define a classical history measure because they are conserved along the 
classical trajectories (recall (13.16bp .) 



IV. HOMOGENEOUS ISOTROPIC MINI-SUPERSPACE MODELS 
A. Euclidean Action and Equations for its Extrema 

From now on, we use Planck units where h = c = G = 1. The Euclidean action /[f?, $] is 
a sum of a curvature part Ic and a part /$ for the scalar field $. The general form for the 
curvature action is: 

Ic[g] = — I (fx{gf'^iR - 2A) + (surface terms). (4.1) 

The general form for the matter action for a scalar field moving in a quadratic potential is: 

U[g, $] = ^ / dHaf'^"^^? + ^'$']- (4.2) 

The integrals in these expressions are over the manifold M with one boundary defining the 
NBWF (cf. eq (11. 2p ). With a convenient overall scale, the homogeneous, isotropic metrics 
are defined as in (13. ip . With that normalizing factor the scale factor a(A), nor the lapse 
N{X), nor any of the coordinates carry dimensions'^. 

It proves convenient to introduce dimensionless measures H, 0, and /i of A, $, and m 
respectively as follows: 

= A/3, (4.3a) 



^ The scaling of the metric used here is different from that employed in [12| . as are others in this paper, 
but they prove convenient for simplifying the numerical work. 
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= (47r/3)^/2^, 
12 = 



(4.3b) 
(4.3c) 



The scaling for H was chosen so that the scale factor of a classical inflating universe is 
proportional to exp{Ht) — the usual definition of H. The other scalings were chosen to 
make the action simple. In these variables the Euclidean action takes the following simple 
form: 



37T 

UP 



dXN 



N 



- a + + a? 



N 



2i2 



(4.4) 



where ' denotes d/d\ and the surface terms in (14.11) have been chosen to eliminate second 
derivatives. The center of symmetry SP and the boundary of the manifold M have arbitrarily 
been labeled by coordinates A = and A = 1 respectively. 

Three equations follow from extremizing the action with respect to A^, 0, and a. They 
imply the following equivalent relations: 



— - 1 + a" 



2 J,2 



0, 



(4.5a) 



am \ N 



1 fa'V 



N \N 



2a 



N 



+ a(l + /xV) =0 . 



(4.5b) 



(4.5c) 



These three equations are not independent. The first of them is the Hamiltonian constraint. 
From it, and any of the other two, the third follows. 

From (14.41) we can read off the explicit forms of the factors in the general form of the 
actions (^Ji) and (jS^D- We have K = 3-k/2H^ and 



Gab = diag(-a, a'), V = (l/2)(-a + i2'(f)'). 



(4.6) 



B. Complex Contours for the Action 



The extremizing solutions a(A), 0(A), and A^(A) will generally be complex. Assuming 
they are analytic functions, the integral (14. 4p can be thought of as taken over a real contour 
in the complex A plane between and 1. Following Lyons [13] it is then useful to introduce 
a new complex variable r defined by 

r(A) = r d\'N{\'). (4.7) 
Jo 

The function r(A) defines a contour in the complex r-plane for each lapse function A^(A). 
Conversely for each contour starting at r = 0, (14. 7p defines a multiplier A^(A) = dT{\)/dX. 
The action (14.41) can be rewritten as an integral over the countour C{0,v) in the complex 
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r-plane corresponding to the A^(A) in (14.41) and connecting r = with an endpoint we denote 
by V. Specifically, 



/[a(r), 0(r)] = / dr \-aa^ -a + a' + a'U'' + ii'^cA] (4.8) 



and / denotes df /dr. 

The equations (14.51) also simplify in the new variable, viz: 



d^ + + (^-0^ + j = 0, (4.9a) 

+ 3(a/a)0-/i^0 = 0, (4.9b) 

a + 2a02 + a(l + /i^^^) = . (4.9c) 

These are the equations we will use to calculate the complex extremizing geometries and 
matter field configurations. 

Using these equations the value of the action (14.81) on a solution can be rexpressed as 

mr), 0(r)] = ^ / dra [a' (l + /i^') - l] , (4.10) 

^-n Jc(0,v) 

Two contours that connect the same endpoints in the r-plane give the same value for 
the action provided they can smoothly be distorted into one another. They are different 
representations of the same extremum as far as the semiclassical approximation to the NBWF 
is concerned and we count their contributions only once. Another way of saying this is that 
(14. 7p defines a complex transformation of the coordinates in the formula for the action under 
which it is invariant if the contours can be smoothly distorted into one another. It should 
not, however, be thought of as a transformation of the coordinates on the manifold which 
remain real throughout. 

This suggests that a solution to equations (14.91) should be considered as a pair of complex 
analytic functions a(r) and 0(t). We can evaluate the action with these functions by picking 
any convenient contour in r connecting the center of symmetry to the boundary. We will 
exploit this in what follows. 



C. Lorentzian Equations 

For the semiclassical approximation to the NBWF we will be interested in complex so- 
lutions to equations (14. 9p . But the ensemble of histories to which these solutions supply 
probabilities will be real, Lorentzian metrics of the form 



df = (3/A) -N\X)dX^ + d\x)dnl 



(4.11) 



We will use hats to distinguish Lorentzian quantities that are always real from the complex 
metrics that extremize the Euclidean action. 

Both the Lorentizian action S and the equations locating its extrema can be obtained 



from the complex relations by substituting = ±iN, a = d and 



into (14.41) . Adhering 
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to the usual convention that the kinetic energy term of the matter be positive, the Lorentzian 
action is 

2 



5[a(A),0(A)] = ^ I rfAiV<j-a[^ 



+ a - + 



212 



We quote the consequent Lorentzian equations in terms of dt = NdX, 

2 



dd 
It 



+ l-d^ -d^ 



dt 



+ n 



212 



0, 



(4.12) 



(4.13a) 



(4.13b) 



d'^d 

dt^ + I dt 



a(l + ji 



2 j,2n 



. 



(4.13c) 



The energy density in the scalar field p$ is a useful quantity for analyzing Lorentzian solu- 
tions. For example, if its exceeds the Planck density we can consider the solution classically 
singular. An expresson for it can be derived from the action (14.121) or from the form of the 
constraint equation fl4.13ap . The result is 



3^ 

Stt 



2 J,2 



(4.14) 



D. The Classical Ensemble of a Complex Extremum 

As discussed in Section IIII Bl the semiclassical approximation to the NBWF given by 
a solution to (14. 9 p corresponds to a solution to the Lorentzian equations (I4.13P when the 
classicality condition (I3.13P are satisfied. The Lorentzian solutions obtained this way are 
the integral curves of 

5(fe,x) = -Im[J(6,x)]. (4.15) 

To calculate this ensemble explicitly for these models we proceed as follows: Choose a 
matching surface of constant 6 = 6,. in a region of minisuperspace where the classicality 
condition is satisfied — typically for large values of 6*. The integral curves of S can be 
labeled by the value of x = X* where they intersect this matching surface. The Lorentzian 
and Euclidean momenta there are given by gradients of S and / respectively [cf. (I3.14bp . 
Their explicit forms in terms of scale factor and field can be found from the actions (14.41) . 
and (I4.12p . From (I4.15P we have on the matching surface: 

6 = 6,, = -Im(pb)|b^, (4.16a) 

X = X*, Px = -^'^iPx)\b*- (4.16b) 

These relations show how a complex extremum specifies Cauchy data for solving the 
equations (I4.13P to find the complete Lorentzian history labeled by x*- The value of 
exp[— 2//j(6*, X*)] gives its relative probability. The classical ensemble is generated as x* 
varies across the matching surface. 
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V. COMPLEX SOLUTIONS 



As the first step in determining the probabihties of the ensemble of classical cosmologies 
predicted by the NBWF we begin by evaluating it semiclassically in this section. The 
classical ensemble implied by this approximation will be determined in the following section. 

We find the semiclassical approximation to the NBWF by numerically solving eqs fl4.9p 
for a(r) and (^(r) along a suitable contour in the complex r-plane connecting the South 
Pole r = with an endpoint v = X + iY where a and take real values b and x- The 
(no) boundary conditions fl3.19al) of regularity at the SP mean that the complex value of 
at the origin is the only free parameter there. To reach the prescribed values (6, x) at 
the boundary, we will adjust both this and the endpoint of integration v. This gives four 
real adjustable parameters to meet four real conditions at v. Hence for each b and x there 
is a unique solution. These solutions can be found analytically in the limits when 0(0) is 
very large and very small. We discuss these limiting cases first as they will motivate our 
numerical search procedure. 

When the scalar field is large (but well below the Planck density) the classical dynamics 
should be governed only by the scalar field and its backreaction on the geometry. The 
background cosmological constant should be largely irrelevant. In this regime therefore we 
expect to recover the approximate complex solutions found by Lyons [12] for a scalar field 
model with a quadratic potential and A = . Indeed, following Lyons one can show that for 
10(0)1 S> 1 the complex Einstein equations (14. 9 p admit the approximate 'slow roll' solution 

Mr) ^ 0(0) + z^, a^r) ^ ^^f'''''^'^^^''^'^'- (5-1) 

There is a similar approximate solution (0_(r), a_(r)) found by changing i to —i in (15. ip . 

These solutions are the complex analogs of the standard 'slow roll' inflationary solutions. 
They can be found up to a constant multiplicative normalization of the scale factor by 
neglecting the cosmological constant, spatial curvature, and terms in equations (14.90 and 
solving the resulting simple set of equations. The results are not good approximations 
everywhere in the complex t = x + iy plane. They hold when when y is not so large that 
the slow role assumption breaks down, and only in regions where |a(r)| >> 1 so that the 
spatial curvature is exponentially negligble as was assumed in deriving them. The existence 
and location of such a region depends on the value of Re[0(O)]. When Re[0(O)] > the 
solution (a+(r), 0+(r)) is valid in a region in the y > half-plane, and (a_(r), 0_(r)) holds 
in a region in the lower half-plane. 

When 10(0)1 is very large and r is sufficiently small that any change in is negligible 
these solutions must match the 'no-roll' solution 

4.(r).m. -ir).?^!^^. (5.2) 

This is regular at the origin and valid in regions in both half-planes. Matching with this 
solution determines the multiplicative normalization factor in (15. ip . 

Following Lyons we now make the further approximation that Re[0(O)] ^ Im[0(O)]. 
Then in the solution (15. ip the scalar field is approximately real along a vertical line r = 
— 30/(O)//i + iy where 0/(0) = Im[0(O)]. Eliminating r from the solution (15.10 for a gives 



2/i0ii(O) 
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where 0i?(O) = Re[(/)(0)]. Therefore, if one takes 



one obtains vertical hnes given by 



TT 

T = -— + iy (5.5) 



along which both a and (p are approximately real. It is clear that progressively finer tuning 
of 0/(0) will yield approximately vertical curves of exactly real a{y) and (f){y). Notice also 
that the condition (15. 4p at the SP that fixes a and to be approximately real does not 
depend on the time parameter along the vertical line. 

The complex action of a solution to the equations of motion can be obtained from (14.101) . 
For the solutions (15.11) the main contribution to the real part Ir comes from the integral 
over real r from the SP to X = tt/ (2/i(/)i^(0)). In this regime the solutions are approximately 
given by the no-roll solution (15. 2p with 0(0) ~ 0_r,(O). When is large this yields 

TT 71 

where we are assuming that at the boundary x ~ 0/?(O). This is the action of Euclidean de 
Sitter space with effective cosmological constant 3m^(0/j(O))^. The main contribution to S 
comes from the integration over the vertical line to an endpoint v. It is given by [l2| 

S - (5.7) 

This can be used to verify whether the solutions satisfy the classicality condition (13.131) at 
large scale factor. One has 

and 

(V^)2^-/iV&'-/iV. (5.8b) 

This means |(Vfe/ij)2|/|(VbS)2| ^ 0. More particularly from (ESI) and (EZI) we have 
\{VM'\/\i^xS?\ ^ l/if^xbr and |(VJ,j)2|/|(V^)2| ^ l/(/i6)V. Hence \{VaIr)'\ « 
{{VASyi, provided X) and hence 0i?(O), is sufficiently large (which we assumed from the 
outset). The complex solutions (15.11) therefore tend to solutions of the Lorentzian Hamilton- 
Jacobi equation along the vertical lines where both a and are real, and they do so in the 
"inflationary" slow roll regime where x is still large. The classicality condition is satisfied. 

Small values of the scalar field are another regime for which analytic approximations 
are possible. These can be calculated by perturbation theory which is the subject of the 
Appendix. There we find that in the leading approximation of vanishing scalar field the 
contours where a(r) is real are exactly vertical in the complex r-plane. In the linear ap- 
proximation the tuning of 7 to give a real value of 0(r) along this contour can be carried 
out explicitly with the result given in Figure [181 

These analytic approximations for large and small scalar field do not, of course, give us 
the complete ensemble of classical histories in these models. To find the complete ensemble 
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FIG. 1: The complex solutions that provide the steepest descents approximation to the NBWF are 
found by integrating the field equations along a broken contour Cb{X) in the complex r-plane. In 
order for the solutions to behave classically at late times one ought to tune the tangent 7 of the 
phase of (p at the South Pole and the turning point X of the contour. We show the tangent 7 (left) 
and the turning point X (right) here as a function of the absolute value (pQ of (/> at the South Pole. 
There is a qualitative difference between < 3/2 models on the one hand where 7 remains finite 
for all (po (dotted curve), and /x > 3/2 models on the other hand where 7 diverges at a critical value 
0g (remaining curves). In the latter case there is no combination {X,'y) for which the classicality 
conditions at large scale factor hold when 4>o < c^qI the ensemble of possible classical histories is 
restricted to a bounded surface in phase space. The right panel shows the critical value (pQ increases 
slightly with for fixed m, and tends to 1.27 as A ^ 0, independently of the value of m^. From 
top to bottom, the different curves show 7 and X for fi = 3/4, 33/20, 9/4, 3 and (in the left panel) 
for a scalar field model with = .05 and A = 0. 



we now solve eqs (14. 9 p numerically in a systematic manner. Guided by the analytic solutions 
(15. ip we begin by taking the scalar field at the origin to be large and approximately real. 
Define (po and 9 by 

<P{0) = \m\e''^<Poe''. (5.9) 

Then with (po ^ 1 and 6 small, we integrate (14. 9 P along a broken contour Cb{X) that runs 
along the real axis to a point X, and then up the imaginary y-axis. We are able to adjust 
both the turning point X and the phase angle 6 so that a and (p tend to real functions b{y) 
and x{y) along the vertical line given hj t = X + iy in the complex r-plane. 

We find that for all fi and for each large (po there exists a unique combination X and 6 for 
which the fields become real at large y. Furthermore for these pairs {X, 6), the ratio of the 
gradients of the real to the imaginary part of the action in different directions all tend to 
zero at large y. Hence the classicality conditions (13.131) hold for this set of solutions, which 
therefore specifies a one-parameter set of classical histories as described in Section IIIII It 
will prove convenient to label the members of this set by (po- 

However, when we decrease (po to (po ~ C>{1) there is an important qualitative difference 
between /i < 3/2 models and models with fi > 3/2. We illustrate this in Figure [1] where we 
plot the values for the turning point X of the contour Cb{X) and for the tangent 7 = tan^ 
of (p at the SP, both as a function of 0o- The values given there have been fine-tuned so that 
the complex solutions behave classically at large scale factor. One sees 7 remains finite for 
all 00 when ^ < 3/2. By contrast in all /i > 3/2 models — and for quadratic scalar field 
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FIG. 2: The real and imaginary part of the scalar field (p for a typical complex solution that 
provides the semiclassical approximation to the NBWF. This solution has /x = 3/4 and (j)Q = 2. 
It is shown along a broken contour Cb{X) in the complex r = {x,y) plane that runs first along 
the X-axis from the South Pole at x = to a value x = X, and then vertically in the y-direction. 
The turning point X is the largest value of x plotted in the left hand two figures. It and the 
imaginary part of (/) at the SP are determined by the requirement that the imaginary part of the 
action becomes constant with increasing y (cf. Figure [T]) . This is necessary for classicality at late 
times and it implies that the imaginary part of (p decays rapidly to zero with increasing y. 



potentials with A = — we find 7 diverges as 0o decreases to a critical value 0q ~ 0{1). 

The classicality conditions (13.131) do not impose a constraint on 00 when fi < 3/2. For 
/i = 3/4, as 00 ~^ 0, it is clear from Figure [1] that 7 — > —.32 and X 7t/2. These limiting 
values agree with the predictions of the perturbation theory for small values of around 
empty de Sitter space, as discussed in Appendix A, FigfTSl 

The behavior of the scalar field 0(r) and the scale factor a(r) along Cb{X) are shown in 
Figures [2] and [3] for a typical complex solution that provides the semiclassical approximation 
to the NBWF. The turning point X is determined by the requirement that the imaginary 
part of the action becomes constant with increasing y (cf Figure [T]). This is necessary for 
classicality at late times and it implies that the imaginary part of and a decay rapidly to 
zero with increasing y. This decay can be exhibited analytically in perturbation theory, see 
Figure [161 

The critical value 0o that is present in all yU > 3/2 models increases slightly with /i, for 
fixed m, and tends to approximately 1.27 when A 0. This limiting value is the critical 
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FIG. 3: The real and imaginary part of the scale factor for the same complex solution as in Figure 
[21 along the same broken contour Cb{X) in the complex r-plane. The imaginary part of a(r) 
rapidly decays to zero along the y-axis as a consequence of the classicality conditions, whereas the 
real part grows exponentially for some time. 

value in standard scalar field models with quadratic potentials and vanishing cosmological 
constant (see also We illustrate this in Fig[T](left) where the solid (black) curve shows 
7 in a A = model with = .05. One has (pQ = 1.27 in this model, and this is independent 
of the mass of the scalar field. 

It also follows from the convergence of the curves in Fig [1] that at large (po, X ~ 1/0^(0) 
and 7 ~ —1/(0^^(0))^ independently of n. This is in agreement with the behavior fl5.4p 
and (15.51) for the analytic solutions (15. ip . Hence, in this regime, only a small complex part 
at the SP is required to reach real b and x l^ite times. The metric representation of 
the complex geometries along the broken contours thus resemble 'fuzzy instantons', with an 
approximately Euclidean section smoothly joined onto an approximately Lorentzian section. 

Most importantly, we find numerically that in the leading semiclassical approximation 
there are no solutions that obey the classicality conditions (I3.13P at large scale factor, other 
than those shown in Fig [H The one-parameter set of solutions given there completely 
determines the ensemble of classical homogeneous and isotropic cosmologies predicted by 
the NBWF. Histories other than these have zero probability in this approximaton. This 
means that whereas the NBWF slices through the whole of phase space for ^ < 3/2, in 
^ > 3/2 models the ensemble of possible classical histories is restricted to a surface in phase 
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FIG. 4: Le/it panel: The ratio of the gradient squared of the real to the imaginary part of the action 
plotted along the y-axis, for the complex solution that behaves classically at large scale factor, with 
/i = 3/4 and = 2. For y < 2 the ratio still significantly deviates from zero, not because Ir 
varies in the y-direction (as can be seen in the right panel) but because the gradient of Ir in the 
X-direction is not small. Right panel: The real part of the action of this complex solution rapidly 
stabilizes along the y-axis. 



space with boundary. 

To demonstrate that the solutions given in Fig [T] satisfy the classicality conditions fl3.13p 
we show in Figure ll]( left) the ratio | V/rP/| VSp for a typical solution for /i = 3/4, along the 
vertical line t = X + iy where a and tend to real functions b{y) and xiv)- As discussed in 
Section [mB] it is necessary for classicality that this ratio tend to zero. We have also verified 
at a selection of points that the ratios of the projections of the gradients, both in the X 
and Y directions, similarly tend to zero at large y. We conclude therefore that this set of 
solutions behaves classically at large scale factor. 

The real part Ir of the action rapidly tends to a constant along the vertical lines where 
the classicality condition holds. This is illustrated in Figure H] (right) for a typical solution 
with yU = 3/4. This means that along these lines the complex solutions become integral 
curves of S. They are the classical Lorentzian histories. The value of the real part of the 
action provides the relative probability of the different classical histories predicted by the 
NBWF. 

The asymptotic value of Ir is shown in Figure [5] as a function of (po and for /i = 3/4. The 
small 00 solutions can be interpreted as classical perturbations of de Sitter space, and can 
be understood analytically in perturbation theory. The upper curve in Figure [5] shows the 
values for Ir predicted by the perturbation theory for small around empty de Sitter space, 
discussed in Appendix A. One sees that this provides a good approximation for 0o < 3/4. 

As mentioned earlier, when ^ > 3/2 there is a critical 0q ~ 0{1) at which 7 diverges. 
Furthermore there are no solutions that obey the classicality condition fl3.13p at large scale 
factor for 0o < 00- In particular, although it is still possible to tune the angle 6 such that a 
and are simultaneously real at some endpoint v = X + iY in the complex r-plane, we find 
the ratio of the projection in the X-direction of the gradients of the real to the imaginary 
parts of the action is always at least of 0{1). 

An analytic analysis of the complex solutions in perturbation theory for small values of 
around the empty de Sitter solution supports this conclusion: In Appendix A it is shown 
that for /i > 3/2 the real part of the action of regular complex solutions of the scalar field 
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FIG. 5: The asymptotic value of the real part of the action of the complex solutions that behave 
classically at large scale factor plotted as a function of (pQ and for fi = 3/4. This determines the 
relative probabilities predicted by the NBWF for the corresponding classical Lorentzian histories. 
The upper curve shows the prediction of the perturbation theory for small (j) around the empty de 
Sitter space with cosmological constant A. 

perturbation equation does not approach a constant along the vertical integral curves of the 
putative classical Lorentzian solutions. Instead it oscillates along these curves, as shown in 
Fig [231 The curves along which 1^ is constant oscillate around some mean value X at large y 
and are shown in Figure EH Along these curves however the ratio of the gradients projected 
in the x direction does not become small at large Y, as shown in Figure [221 Hence, at least 
in the semi classical approximation, all classical histories for small scalar fields in all /i > 3/2 
models have zero probability in the NBWF. Even in perturbation theory the classicality 
condition is non-trivial. 

For large values of the scalar field the background cosmological constant is largely irrele- 
vant in the early universe. In this regime the complex extremizing solutions are qualitatively 
similar in all models we have considered. The asymptotic value of the real part of the action 
of the complex solutions that imply classical behavior at late times is shown in Figure [5] 
as a function of 0o in three /i > 3/2 models (left) and in a A = model with a quadratic 
potential (right). The action tends to a finite value at the lower bound (pQ that arises from 
the classicality condition. At large 0o it goes to zero as ~ 1/0^(0), in agreement with the 
behavior (15. 6p that follows from the slow roll approximation. One sees that without further 
constraints the NBWF universally favors histories with (po near the lower bound (pQ. 

VI. CLASSICAL HISTORIES 

The complex 'fuzzy instantons' that extremize the Euclidean path integral defining the 
NBWF were discussed in the previous section. They provide the probabilities for the en- 
semble of real Lorentzian classical histories which are the subject of this section. 
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FIG. 6: The asymptotic value of the real part of the action of the complex solutions that behave 
classically at late times plotted as a function of (^Qj in three /i > 3/2 models (left) and in a scalar 
field model with quadratic potential and A = (right). The values of /i and rv? are as in Figure [TJ 
The action tends to a finite value at the lower bound 0g that arises from the classicality condition, 
and it goes to zero as ~ 1/0^(0) at large </>o. 



Complex extremizng solutions that predict classical behavior obey the no-boundary con- 
dition at the SP and the classicality condition fl3.13p at the boundary at large scale factor. 
The values of a and together with their derivatives at the boundary provide Cauchy data 
for the ensemble of classical Lorentzian histories predicted by the NBWF as discussed in 
Section lIVDi In this section we study various properties of the members of this ensemble 
by evolving these Cauchy data backwards and forwards in time using the Lorentzian field 
equations fl4.13p . Combined with the results for the relative probabilities provided by the 
action of the complex solutions (e.g. Figs [5] and [6]) this allows one to predict probabilities for 
several features of our specific universe if it is in the no-boundary state. These include the 
amount of inflation, whether it had an initial bounce or singularity, its future behavior, its 
time asymmetry if bouncing, and its consistency with the standard cosmological model . We 
will continue to label the individual classical histories in the ensemble by 0O) the absolute 
value of the scalar field at the SP of the corresponding complex solution. 



A. Inflation 

For large 0o the complex solutions are well approximated by the analytic form (15. II) . 
Moving upwards along the vertical contours where a and are real and 0o ^ 1 has 
(with y{t) = t) 

m = mt)) ^ MO) - f , a{t) = a{y{t)) ~ e^<^W*, (6.1) 

until the scalar field becomes less than ~ 1/2. This is just like the behavior of Lorentzian 
slow roll inflationary solutions, and it shows the classicality condition implies inflation at 
large 00 . 

This connection is general. Figure [7] shows the trajectories of several numerically 
calculated histories in (/i, 0) variables, where h{t) is the instanteous Hubble constant 
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FIG. 7: The no-boundary wave function predicts that ah histories that behave classicahy at late 
times undergo a period of inflation at early times as shown here by the linear growth of the 
instantaneous Hubble constant /i = a^f/a in five representative classical histories for // = 3 and for 
(j)Q between 1.3 and 4. 

h = [da / dt) / a = d^t/d. Five representative members of the ensemble of classical histo- 
ries for n = 3 and for 0o between 1.3 and 4 are shown. When we follow the histories back 
in time to higher values of h and 0, they all lie within a very narrow band around h = fi(f). 
This is characteristic of Lorentzian slow roll inflationary solutions. Furthermore, since the 
numerical analysis shows that there are no solutions other than those given in Fig [1], and 
represented in Fig [71 we conclude that the NBWF predicts that a classical homogeneous 
and isotropic universe must have an early inflationary state. The NBWF and classicality at 
late times imply inflation at early times. This conclusion holds for all values of fi. 

Although the classicality condition implies inflation for all values of /i the drivers of 
inflation are different for different values. For /z < 3/2 and small 0o inflation is always 
driven by the background cosmo logical constant. In all other models however, and for small 
fi at large (pQ, inflation in the early universe is driven by the scalar field potential energy. By 
this we mean specifically that d^u > when > .5, which we find is the minimum value of 
required for inflation to occur in A = models with quadratic potentials and small m^. 

To get a quantitative measure of the amount of inflation predicted we calculated the 
number of efoldings N = J hdt of scalar field driven inflation over the range of time where 
d^tt > and > .5 for the members of the ensemble of Lorentzian histories predicted by the 
complex solutions found in Section 5. The results are summarized in Figure [H For /i > 3/2 
the lower bound on 0o arising from the classicality condition implies the number of efoldings 
is always greater than one^ . It follows from Figure [HI combined with the information on the 
relative probability of the histories given in Fig [5] and Fig [6l that by itself the no-boundary 



^ For initially singular solutions inflation generally does not begin immediately at the singularity, see e.g. 
Figure [Til 
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FIG. 8: Left panel: The number of efoldings N of inflation driven by the scalar field (as opposed to 
the cosmological constant of the background) in the classical histories predicted by the NBWF, for 
five different models. From top to bottom, the different curves correspond to /x = 3/4, 3, 9/4, 33/20 
and A = 0, = .05. Right panel: Detail of the left panel showing the regime around the critical 
value (f)Q in the /.i > 3/2 and pure scalar field models. The lower bound on (f)Q that arises from 
classicality implies a lower bound on the number of efoldings. 

wave function favors a small number of efoldings on a history by history relative probability 
basis. The answers to more physical questions involving probabilities conditioned on the 
data in our past light cone are obtained from the no-boundary probabilities by summing 
them over those for classical spacetimes that contain our data at least once, and over the 
possible locations of our light cone in them. This sum can significantly change the no- 
boundary predictions based on the wave function alone 0, [2l[ . We return to this point in 
Section Ivni 

B. Bounces and Initial Singularities 

For /i > 3/2 and — ^ allowed classical histories of the universe are singular at 

an initial time tg. Near the singularity both the potential and the curvature are unimportant 
in the Einstein equations fl4.13p . and one has a{t) ~ {t — tsY^^ ~ ln(t — t^). But for 

00 > 00 histories bounce at a finite radius at in the past. The critical value 0o at which 
there is a transition from initially singular to bouncing is determined entirely by the Einstein 
equations, and therefore independent of H. A bounce at a finite radius in the past is possible 
despite the singularity theorems because a scalar field and the cosmological constant violate 
the strong energy condition. Even though such non-singular classical solutions form only 
a small subset of all scalar field gravity solutions they have significant probability in the 
no-boundary state. Near a bounce the universe approaches a de Sitter state with radius 
~ {Hfi(f)b)~^ where 0f, is the value of the scalar field at the bounce. For sufficiently large 0o, 
0b ~ 00; as shown in Fig [9] (right panel). The scale factor at the bounce versus 0o is shown 
in Fig [9] (left panel), which clearly reveals the transition from bouncing solutions to initially 
singular ones. The critical value 0o itself is shown as a function of fi in Fig [101 One sees 0o 
slightly increases with /i asymptoting to ~ 1.54 as A — > 0, for fixed m. As discussed for 0g 
above, this limit corresponds to the critical value that separates the bouncing from singular 
histories in pure scalar field models with quadratic potentials for any nonzero mass < 1. 

When we evolve the data provided by the complex solutions on the matching surface 
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FIG. 9: The scale factor en, (left) and the scalar field 05 (right) at the bounce of the classical 
histories predicted by the NBWF. The values of ^ and A are as in Figure El When /j, < 3/2 the 
histories always bounce at a minimum radius in the past. By contrast for n > 3/2 there is a 
transition from bouncing to initially singular at a critical value (pQ ~ 1.5. Above that (pb ~ <Ao- 
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FIG. 10: A 'phase diagram' that summarizes some of the early universe properties of the Lorentzian 
histories predicted by the NBWF. For /x < 3/2 there is one classical history associated with each 
value of (j)o and the universe bounces in the past for all ranges of (pQ. For fi > 3/2, however, 
there are no classical histories for (po < (pQ (bottom curve). Between cpQ and (pQ (top curve) the 
Lorentzian histories have an initial singularity. Finally, for (pQ > (pQ the Lorentzian solutions bounce 
at a nonzero minimum radius in the past. The limiting values of cp^ and cp^ when A ^ 0, for fixed 
m, are 1.27 and 1.54 respectively. 

backwards in time we find that, for < 3/2, all the histories bounce at a minimum radius 
in the past. Hence, in this regime, the classicality conditions at late times select a set of 
histories in the NBWF where either the potential energy of the scalar field or the background 
cosmological constant dominate the evolution at early times. 

Figure [11] shows the ratio of the scalar field energy density (14.141) over the vacuum energy 
density at the bounce, again as a function of 0o- Since 06 ~ 0o for most 0o; (Figure [H]) this 
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FIG. 11: The ratio of the energy density pf, in the scalar field at the bounce over the vacuum 
energy density py^ for the allowed Lorentzian histories for [i = 3/4 (left) and for (right, from top 
to bottom) // = 3, 9/4 and 33/20. When /x > 3/2 there is a minimum matter density needed at 
early times for the history to exhibit classical behavior at late times 
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FIG. 12: A 'phase diagram' that combines some of the early and late time properties of the 
Lorentzian histories predicted by the NBWF. Below the solid curve the homogeneous isotropic 
classical universes predicted by the NBWF recollapse to a Big Crunch, whereas above the curve 
the universes continue to expand forever. 



generally grows quadratically with {/)o. When /i < 3/2 we find classical histories for all ranges 
of densities, whereas for yU > 3/2 there is a minimum matter density needed for the history 
to exhibit classical behavior at late times (except for the vacuum de Sitter solution). For 

realistic values of A and fi, therefore, a nearly empty de Sitter solution has zero probability 
in the semiclassical approximation to the NBWF. 
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FIG. 13: The reality of the NBWF implies the ensemble of ahowed classical histories is time- 
symmetric. The individual classical bouncing histories, however, are generally not time-symmetric 
about the bounce. A natural measure of the amount of time-asymmetry of the individual histories 
is provided by the quantity rj = [4>^t)b/4>b- The left panel shows 77, as a function of for /i = 3/4. 
In the right panel we plot r/ for (from top to bottom) /i = 33/20, 9/4 and 3. One sees that for 
/i > 3/2, r/ diverges when (pQ 0g where (/>q separates singular from bouncing solutions. 



C. Eternal Expansion and Final Singularities 

Next we explore the late time properties of the classical histories. Evolving the data 
provided by the complex solutions on the matching surface forward in time, the universe 
either expands forever or recoUapses again to a singularity^. We find that when /i < 18 the 
histories expand forever for all ranges of 0o that admit classical histories. In these universes, 
when the scalar field rolls down the potential, the cosmological constant of the background 
takes over to drive the exponential expansion. When /z > 18, however, there is critical value 
0Q, and for (pQ < 0q the universe recollapses. We plot 0g as a function of /i in Fig [12], which 
shows that this slowly increases with fi. 



D. Time Asymmetry 

Bouncing classical histories are generally time-asymmetric about the bounce. This can 
be seen in perturbation theory (e.g. Fig [H]) and a particular non-perturbative case is the 
fourth example in Figure O A natural measure of the time asymmetry of histories at the 
bounce is given by 

V = {kt)b/k- (6.2) 

We plot ?7, as a function of 0o in Figure [131 fo^' = 3/4 (left) and for several values of 
/i > 3/2 (right). In the former model, the limiting value 77 ^ .11 for 0o ^ agrees with the 
prediction of the perturbation theory for small 0, which we obtain in Appendix A. One sees 
that, for /i > 3/2, r] diverges when 0o — ^ 0o — boundary between singular and bouncing 
histories. For large 0o- ''7 — ^ in all models. 



We did not find any Lorentzian solutions with multiple bounces 
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However, our numerical search 



procedure identifies discrete solutions and may miss such highly fine tuned examples. 
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At the current level of our analysis in which we restrict attention to homogeneous isotropic 
minisuperspace models it is not clear whether the time-asymmetry of the bouncing classical 
histories in the NBWF has any physical (observable) effects. But one might expect ob- 
servable signatures of the time asymmetry to show up in the spectrum of inhomogeneous 
perturbations. We intend to calculate these in future work. We emphasize also that although 
individual classical bouncing histories are not generally time-symmetric about the bounce, 
the reality of the NBWF implies the ensemble of allowed classical histories is time-symmetric. 
For every history in this ensemble, its time reversed is also a member. 



E. Cosmological Models 

We next turn to particular kinds of Lorentian histories. A gallery of qualitatively different 
classical histories for one value of 0o is exhibited in Figure [TH A class of histories that are 
particularly interesting from an observational point of view are represented by points that 
lie just above the curve of (pQ vs (po over a range of /i in Figure [121 These turn out to 
correspond to universes that undergo an early period of inflation that is followed by an era 
of oscillating scalar field that will lead to matter generation and domination. Eventually 
the cosmological constant takes over to drive a second phase of exponential expansion that 
lasts forever. The NBWF, therefore, appears to be consistent with the standard picture 
of inflationary cosmology for our universe in which a scalar field rolls down from high up 
the potential and subsequently oscillates around the minimum losing its energy into created 
particles. An example of a history of this kind is given in Figure [Ml (2nd row). 

Further increasing fi for fixed 0o yields a qualitatively different universe. For /i > 98 
the 00 = 1-32 histories lie below the solid curve in Fig [121 These universes have an initial 
singularity and recollapse again to a Big Crunch. An example of a universe of this kind is 
given in Fig[TH 1st row. The third example in Fig[T^ shows the (initially singular) 0o = 1-32 
history for /i = 9/4, where the cosmological constant immediately takes over to drive the 
expansion when the scalar field has rolled down its potential. Inflation never really ends in 
this universe. 



VII. VOLUME WEIGHTING 

The NBWF gives the probabilities of entire classical histories. But we are interested in 
probabilities that refer to our data, which are limited to a part of our past light cone. Among 
these are the top-down probabilities for our past conditioned on (a subset of) our present 
data [24]. Hawking [2l[ and the current authors [3| have argued that in homogeneous models 
these are obtained by multiplying the NBWF probabilities for classical histories by a factor 
exp(3A^) proportional to the volume of the hypersurface on which our data approximately 
lie^. This multiplication can be understood as resulting from a sum over the probabilities 
for classical spacetimes that contain our data at least once, and over the possible locations 
of our light cone in them . In a large universe there are more places for our data to be. 

In order for the volume weighted probabilities to be physically meaningful as probabilities 
relevant for what we observe the universe must obviously last to the present age of 14 Gyr. 



^ Anthropic reasoning has also been used as an argument to include a volume factor 25|. 
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FIG. 14: A gallery of classical Lorentzian histories. The scale factor d{t) (left) and the scalar 
field (p(t) (right) in the classical histories labeled by (po = 1.32, for four different values of ji. The 
value of fi decreases from top to bottom taking the values fi = 100,96,9/4,33/20. If is fixed 
these correspond to increasing A. These are four qualitatively different cosmologies, ranging from 
initially singular histories that recollapse again (top) to eternally expanding universes that bounce 
in the past (bottom). At intermediate values of fi (2nd row) the NBWF is consistent with the 
standard picture of inflationary cosmology, consisting of a short period of inflation as the scalar 
field rolls down from high up the potential followed by an era of oscillation representing particle 
creation and ensuing matter domination. Eventually the cosmological constant takes over to drive 
a second (future-eternal) phase of exponential expansion. 
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FIG. 15: To account for the different possible locations in the universe of the Hubble volume that 
contains our data, one ought to multiply the relative probabihties for classical histories coming 
from the NBW (left) by a volume factor, to obtain the probability (right) for what we observe in 
our past light cone. The resulting volume-weighted probability distribution favors a large number 
of efoldings when the ensemble is restricted to universes that last sufficiently long. 



This further restricts the ensemble of histories, requiring (f)Q to be larger than a critical value 
0Q ~ 2 or, equivalently, > 5. 

Figure fTSl shows the qualitative effect in m^0^ models of multiplying the relative probabil- 
ities for classical histories exp(— 2/r) coming from the NBWF by a volume factor exp(3A^). 
Volume weighting clearly enhances the probability for a large number of efoldings. An im- 
portant feature of the volume weighted probability distribution is that there is a wide region 
where the probability is strongly increasing with A^. Indeed when one considers the prob- 
ability distribution ~ exp(3A^ — 21 n), as a function of the value = 0i at which inflation 
starts, the gradient of this probability distribution is positive provided 

> \VJ^- (7.1) 

For quadratic potentials this condition is satisfied well below the Planck density. 

For a realistic value of m Figure [15] shows qualitatively that the two constraints of clas- 
sicality and minumum age yield a restricted ensemble of histories whose volume weighted 
probabilities slightly favor a large number of efoldings. This can be understood analytically. 
Since —Ir ~ l/(m0ij(O))^ ~ l/(m0o)^ for the slow roll solutions predicted by the NBWF, 
and since (pi ~ 4>r{0) ~ 0o, the volume factor exp(3A^) is comparable to the no-boundary 
weight exp(— 2Jij) for 0j = 1/m, i.e. for solutions that start inflating near the Planck den- 
sity. Hence the volume weighted probability distribution is peaked both at low 0o and for 
solutions that start inflating near the Planck density. The latter peak slightly dominates 
when the constraint that the universe lasts ~ 14 Gyr is taken in account (solid curve). 

We expect the effect of the volume factor on the probability distribution to be much 
more dramatic in the context of a landscape potential [s^. Indeed it appears likely that in 
some regions of a landscape potential and in particular around broad saddle-points of V, 
the volume factor more than compensates for the reduction in amplitude due to the higher 
value of the potential. This would lead to the prediction that in a landscape potential, the 
most probable universe consistent with our data had a large number of efoldings and began 
in an unstable de Sitter like state near a broad saddle-point of the potential. Because the 
dominant saddle-points are well below the Planck density we furthermore expect that the 
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most probable histories lie entirely in the semi-classical regime 2J, |26l . 



VIII. ARROWS OF TIME 

Suppose that our classical universe bounced at an early time at a radius well above the 
Planck length. At no time in its history were there large quantum fluctuations in the geome- 
try of spacetime. Could events, structures, and processes before the bounce have influenced 
events, structures and processes today? Could we receive information from intelligent aliens 
living before the bounce encoded in gravitational waves, neutrinos, or boxes made of some 
durable form of matter not yet discovered by us? 

The overwhelming observational evidence for an early hot period in the universe suggests 
that most information from before the bounce could not get through to us in any accessible 
form. Matter was in thermal equilibrium at least at temperatures high enough to dissociate 
nuclei and information encoded at lower energy scale phenomena would be wiped out. But 
gravitational waves whose coupling to matter is the same as that governing the expansion 
may not participate in this equilibrium. 

But even if information could propagate from one side of the bounce to another we have 
to consider the thermodynamic arrow of time to discuss whether events on one side could 
influence events on the other. Causation is generally possible only in that direction. That, 
for instance, is why we remember past events but not future ones. 

In the trenchant analysis of the arrow of time by Hawking, LaFlamme and Lyons [l^ the 
thermodynamic arrow of time is taken to coincide with the time direction in which fluctua- 
tions away from homogeneity and isotropy grow. Small fluctuations grow under the action 
of gravitational attraction into large inhomogeneities. That is order into disorder. Hawking, 
LaFlamme and Lyons examine the evolution of fluctuations in the extremizing solutions that 
provide the semiclassical approximation to the NBWF. They show that regularity conditions 
at the South Pole imply that the fluctuations in the extremizing solutions are small there 
and therefore increase away from the South Pole because they have nowhere to go but up. 

The Lorentzian histories predicted by the NBWF are not the same as the extremizing 
solutions, but they are closely connected. In the homogeneous, isotropic case for example 
the curves of constant real part of the action are Lorentzian trajectories in minisuperspace 
when the classicality condition fl3.13p holds. We have not yet calculated the fluctuations to 
these models. But when we do it seems reasonable to suppose that regularity at the South 
Pole will imply that the fluctuations are small near the bounce and tend to increase away 
from it for a signiflcant time. 

Assuming this result, the arrow of time in bouncing solutions increases away from the 
bounce^. Put differently it points in opposite directions on opposite sides of the bounce. 
It therefore seems unlikely on general thermodjTiamic grounds that events on the opposite 
side of the bounce could influence events on this side. To do so their influence would have 
to travel backward in time. Unless intelligent aliens flnd some way to send information 
backward in time over billions of years we are as unlikely to flnd any messages from them 
as we are to find ones sent by intelligent aliens in our own future. Can we say then that 
the other side of a bounce is 'real' ? It is just as real as the pocket universes in an eternally 
infiating spacetime which also can neither communicate with us nor influence us. 



This has also been considered by CarroU and Chen in a different context 27 1 
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This situation is in sharp contrast with the causahty in ekpyrotic cosmologies 28|] and in 
the 'pre-big bang' models discussed in [29l| , where one typically starts with an ordered state 
in the infinite past and as the universe evolves, departures from this state grow in time. 
Hence in these models the arrow of time always points forward. 

However, a subset of the class of histories predicted by the NBWF has a singularity in 
the past. We have seen (cf. Fig [T5]) that the time asymmetry becomes infinitely large as 
we approach the regime of initially singular solutions. The NBWF does not tell us whether 
evolution continues past this singularity, and it has in fact not been shown rigorously whether 
this is possible at all in any realistic model^. But if evolution continues past this singularity, 
it is conceivable based on Fig [13] that the arrow of time in these histories will always point 
in the same direction. This subset of histories may therefore represent the pre-big bang 
spacetimes predicted by the NBWF in which the arrow of time always points forward and 
information can propagate from the contracting phase to the expanding regime. 



IX. CONCLUSIONS 

The large scale properties of our specific universe can be summarized in a short list 



of facts [31|: Classical physics applies on coarse-graining scales above the Planck length. 
The universe is expanding from a hot big bang in which light elements were synthesized. 
There was a period of infiation, which led to a fiat universe today. Structure was seeded 
by Gaussian irregularities, which are relics of quantum fiuctuations. The dominant matter 
is cold and dark, and there is dark energy which is dynamically dominant at late times. 
Very roughly this list of features constitutes the standard cosmological model. Quantum 
cosmology seeks to provide a theory of the quantum state of the universe that would predict 
connections between these facts. 

The first item on the list — the wide range of time, place, and scale on which classical 
physics applies — is central to all the others. This quasiclassical realm is such a manifest fea- 
ture of our experience that most treatments of cosmology assume it. But, classical behavior 
is not a general feature of quantum systems. Rather, it emerges only for particular coarse- 
grainings in a restricted class of states. That is especially true for the emergence of classical 
spacetime geometry in a quantum theory of gravity^°. Any viable theory of the quantum 
state of our universe must predict classical spacetime over the whole of its visible part from 
the Planck epoch to the distant future. Broadly speaking, this paper has mainly focussed 
on two issues connected with the emergence of a classical cosmological spacetime from the 
'no-boundary' theory of its quantum state: (a) What is the ensemble of classical histories 
predicted by the NBWF and what are their probabilities? (b) What are the implications of 
the classicality condition for the standard model of cosmology? In particular, what are the 
important properties of the members of the ensemble of classical histories predicted by the 
NBWF, and what are the resulting probabilities for what we observe in our past light cone? 

We have analyzed these issues in a very simple class of homogeneous, isotropic minisuper- 
space models with a single scalar field moving in a quadratic potential and a cosmological 
constant. Our main results are as follows: 



for recent work on this. 



^ See however 

Eternal inflation is sometimes said to vitiate the dependence of the present universe on the details of its 



initial quantum state. But those statements typically assume that spacetime geometry is classical. 
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Classical Prediction: Generalized quantum mechanics is a clear framework for the predic- 
tion of classical behavior from the NBWF [1, [l5|. Probabilities are predicted for an ensemble 
of four-dimensional classical histories of geometry and matter field. The complex 'fuzzy in- 
stanton' metrics that extremize the sum-over-histories defining the NBWF are distinct from 
the real Lorentizan classical metrics for which they provide the probabilities. 

The no-boundary measure of the universe: The probabilities for histories in the NBWF 
classical ensemble define a measure on classical phase space. The NBWF measure is concen- 
trated on a surface in phase space which in realistic models has a boundary arising from the 
classicality condition. It is this concentration to a bounded surface in phase space that gives 
the NBWF predictive power. More specifically, for given fi the NBWF generally singles out 
at most a one parameter subfamily from the two parameter family of classical, Lorentzian 
homogeneous, isotropic solutions. For fi < 3/2 we found classical histories for all ranges 
of possible matter content. But for the more realistic case of > 3/2 we found a certain 
amount of matter (0o > 1-27) is necessary for classical behavior if there is any matter at 
all. For yU > 3/2 a nearly empty, almost deSitter space is not the most probable Lorentzian 
history. A significant amount of matter is required for classical histories. 

Inflation and Classicality: All allowed histories that behave classically at late times inflate 
at early times near the bounce or the initial singularity. For /i < 3/2 and small 0o the 
cosmological constant drives the inflation. For /i > 3/2 the required matter is the driver. 
The NBWF and classicality imply inflation. This result illustrates the predictive power of 
the NBWF. Indeed, using a measure extending over all of phase space motivated by classical 
dynamics Gibbons and Turok found a negligible probability for inflation [lo| . 

Number of efoldings: As Fig O and Fig El combined with Fig [HI show, the NBWF on 
its own favors a small number of efoldings on a history by history relative probability basis. 
However, we can ask the more physical, top-down, question of what is the probability of 
inflation in our past conditioned on our limited present data in a Hubble volume. Then the 
probability for a long period of inflation is enhanced as discussed in Section IVIII Roughly 
inflation leads to a larger universe with more possible locations for our Hubble volume. 
Requiring that the universe lasts to the age of 14 Gyr inferred from observation enhances 
the probability for a long period of inflation further. 

Bounces and Initial Singularities: Some histories of the NBWF classical ensemble bounce 
at a minimum radius and some are initially singular. The diagram in Fig [TU] shows the range 
of parameters corresponding to each. On a history by history relative basis Fig O and Fig O 
show that the NBWF prefers singular beginnings. But we can again ask the more physical, 
top-down question of what is the most probable origin of the universe given our limited 
present data in a Hubble volume. Then, as discusssed in j3| and in Section IVIII here, the 
most probable origin may be (depending on the model) a bouncing universe in which the 
universe was always in the semi classical regime. 

Future- eternal expansion and Final Singularities: The NBWF predicts probabilities for 
classical histories and therefore for their long term fate just as much as for their origins. 
Recollapse to a singularity and future-eternal inflation are the two possible futures for ho- 
mogeneous models. Figure [12] shows the range of parameters fi and 0o that correspond to 
each. Recollapse is possible only for large /i (small A) and for universes that have an initial 
big bang singularity as well as a final singular big crunch. 

Singularity Resolution: Even for classical histories that are singular at early times the 
NBWF unambiguously predicts probabilities for late time observables such as CMB fluctu- 
ations. That is because it predicts probabilities for histories rather than their initial data. 
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The NBWF therefore resolves the big bang singularity, in the sense that it is no longer an 
obstruction to prediction. 

Time Asymmetry and the Arrow of Time: Figure [13] shows that individual Lorentzian 
histories are generally time asymmetric although the ensemble of histories is time-symmetric 
on general grounds. For large 0o this asymmetry is small. The restriction to homogeneous 
models does not permit a conclusive discussion of the thermodynamic arrow of time. How- 
ever, one possibility is that it points away from a bounce on either side. Causality in this set 
of histories would be very different from causality in ekpyrotic and pre-big bang cosmologies, 
where the arrow of time always points in the same direction. 

There is much to be done to extend these models to more realistic ones and to back 
up various theoretical assumptions that have been made. Two extensions are of particular 
importance: First, relaxing the restriction to homogeneous and isotropic models would allow 
consideration of the evolution of quantum fluctuations whose effects could be detectable in 
the CMB. Further, bubble nucleation, and the arrow of time could be discussed. Second, as 
we suggested in |3| in more realistic landscape potentials the classicality condition can act 
as a vacuum selection principle resulting in top-down probabilities that favor a bouncing 
universe that had a long period inflation and was always in the semi-classical regime. 

On a technical level, going beyond the lowest semiclassical approximation that has been 
used here could yield more satisfactory probabilities and facilitate comparison with other 
measures such as the classical one developed in [lo|. A study of physically realistic coarse- 
grainings of spacetime geometry and the decoherence of sets of alternative histories defined 
by them would help back up a number of assumptions that we have made. 

These opportunities for extension, however, should not obscure the fact that our results 
in the simple models of this paper already demonstrate that the NBWF and classicality 
condition can play a central role in understanding what we observe of our quantum universe. 
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APPENDIX A: PERTURBATION THEORY 

When the scalar field is small it is a perturbation on the model with only a cosmological 
constant. This is not a physically interesting case since we do not live in a nearly empty 
de Sitter space. But, it is a case where the entire discussion can be carried out essentially 
analytically so as to provide a guide for the detailed numerical calculations for the physically 
interesting non-perturbative situations. 
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This appendix discusses the first two orders of perturbation theory. To avoid a debauche 
d'indices we will use a, a for the leading order Euclidean and Lorentzian scale factors driven 
only by a cosmological constant, and 0, for the perturbations in the scalar field. 



1. No-Boundary Semiclassical Solutions 

We first calculate the complex solutions to the equations (14. 9 p that provide the semiclas- 
sical approximation to the no-boundary wave function at given real values of {b,x)- These 
are solutions which are regular at the South Pole and match the given values at the other 
endpoint. 

No Scalar Field: When there is no scalar field our only concern is the geometry. There 
is only one solution of fl4.9ap that is regular at the South Pole r = and that is 



As discussed in Section IIVBI , finding a regular solution for a given value of b means finding 
a contour in the r-plane connecting the origin to a point t = v = X + iY where a{v) is 
real and equal to b. The scale factor a(r) = sin(r) is real along the curves y = and along 
X = ±tt/2,x = ±37r/2, ■ ■ ■ . If 6 < 1 there is a solution with v on the real axis. If 6 > 1 
there are two candidate solutions corresponding to complex conjugate values of v along the 
constant x curves where a{x + iy) is real with y = ± cosh~^(6). We will argue in a moment 
that only X = tt/2 corresponds to a solution on the no-boundary manifold. 

If 6 < 1 the contour between r = and t = v can be chosen to lie on the real axis. Then 
the metric is real, Euclidean and corresponds to part of the Euclidean 4-sphere. For 6 > 1 
consider the solution with X = 7i/2 and Y = cosh~^(6). The connecting contour can be 
taken to run along the real axis to X = 7r/2 and then up the |/-axis to Y. This corresponds 
to the geometry of half a unit radius round Euclidean four-sphere joined smoothly across a 
surface of vanishing extrinsic curvature to half of a Lorentzian de Sitter space starting at 
the bounce. This is the well known no-boundary instanton j20| nucleating de Sitter space, 
and we will call this the NBI contour. 

Proceeding along a contour of real r to X = 37r/2 and then upwards to Y gives the same 
geometry as with X = tt/2 but with an additional Euclidean four-sphere attached at the 
South Pole. This geometry is not strictly regular on the no-boundary manifold. We therefore 
exclude it and all candidate solutions with larger values of X. Solutions with X = —n/2 
are the same as those with X = -|-7r/2. 

The solution with X = n/2 and Y = — cosh^^(6) has the opposite sign of the imaginary 
part of the action from the one with positive Y. We therefore count it as an independent 
semiclassical solution. The two solutions make complex conjugate contributions to the wave 
function ensuring that it is real as discussed in Section IIIIDi These two solutions dominate 
the semiclassical approximation to the no-boundary wave function when there is no scalar 
field. 

Perturbing Scalar Field: The first order scalar field perturbations to the empty extrem- 
izing solutions for the scale factor found above satisfy fl4.9bl) with a(r) given by flAip . The 
boundary conditions are that 0(r) be regular at the South Pole r = (meaning that (p 
vanishes there) and match the given x the boundary. 



We denote by G(r) the (regular) solution to dilgb]) with G(0) = 1 and G{0) = 0. This is 



a(r) = sin(r) = sin(x -|- iy). 



(Al) 



G'(r)=F(a,6,2,(l-cos(r))/2) , 



(A2) 



37 




FIG. 16: The complex extremizing solution G. G{t, Jl) is the complex solution of the linearized 
equation for the scalar field which is regular at the origin and equal to 1 there. All other field 
extrema are multiples of G. Three plots are shown for Ji equalling .25, .50, and .75. Each uses the 
NBI contour in the r-plane that extends horizontally from the origin along the real axis to x = tt/2 
and then vertically in the imaginary (y) direction. Along this particular contour the unperturbed 
metric makes a smooth transition between a Euclidean instanton and a Lorenzian deSitter metric. 
G is real along the real part of the contour but complex along the imaginary component. The real 
parts of G are indicated by solid lines, the imaginary parts by dashed lines. The curves are in fact 
continuous along the contour with appropriate matching conditions for the derivatives reflecting 
the change in direction of the contour at = 7r/2, y = 0. (cf Figure O) 

where F(a, 6, c, z) is the hypergeometric function, /2 = 2/i/3, and 

a=(3/2)(l + v/l-/x2), 

6=(3/2)(l- v^r^) . (A3) 

(Evidently, G depends on as well as r but we will not usually indicate this explicitly.) 
The function G{t) is real on the real axis for x < vr and therefore real analytic, specifically 
G*{t) = [G{t*)]* = G{t). The function G{t) is multi-valued and there is therefore a cut 
which can be taken to extend along the real axis from x = vr to infinity. We will assume 
that we are considering perturbations of the no-boundary instanton solution defined by the 
NBI contour discussed above and that this lies on the first sheet of G{t). The behavior of 
G{t) for several values of Ji is illustrated in Fgures [11] and [T71 

The general regular solution will be a complex number times G{t). This number can be 
found as follows: The value of b determines a point in the complex plane from the zeroth 
order calculation above. For 6 < 1 this is on the real axis at {X = sm-\b) < n/2,Y = 0). 
Thus the required solution is 

0(r; b,x)=X G(r)/G(sin-i(6)), (6 < 1) . (A4a) 

This is not an especially interesting case from the point of semiclassical prediction since the 
imaginary part of the action, S, will vanish. 

The interesting case 6 > 1 is similar. For definiteness, focus on the specific extremizing so- 
lution whose zeroth order approximation is labeled by the point {X = ti/2,Y = + cosli^^(6)). 
The regular solution for the scalar field matching x that value is 

0(r; 6, x) = A G(r)/G(7r/2 + t cosh-i(6)), (6 > 1) . (A4b) 
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This solution is unique once the contour defining the unperturbed solution is fixed. 

For the non-pertubative case discussed in Sections |V] and |VT] analytic solutions are not 
available. Typically we start at the origin with a value for 0(0), integrate along some 
contour in the r-plane, and adjust the endpoint of integration and the complex value of 0(0) 
to reach given real values of b and x- This connection between 0(0; b, x) and (&, x) is given 
in perturbation theory by Eqs flA4p . e.g for 6 > 1, 

(f){0;b,x)=Mb,x)exp[te{b,x)]=x/G{7r/2 + tcosh-\b)). (A5) 

Since x is real we have 

e{b, x) = -Arg[G(7r/2 + i cosh^^(6))] (A6) 
where Arg is the complex phase, and 

Mb,x) = x/\Gi7r/2 + tcosh-\b))\ . (A7) 
Figure [18] shows the perturbative values of 7 = tan(6') for various values of fi. 



2. The Classical Ensemble 

We next turn to the perturbative construction of the ensemble of classical Lorentzian his- 
tories predicted by the no-boundary wave function and the evaluation of their probabilities. 
The prescription for this is described in Section IIVDI and is straightforward to implement 
explicitly in perturbation theory. For each extremizing solution, we chose a matching sur- 
face 6 = 6^, in minisuperspace and at each point along it labeled by x = X* we evaluate 
initial data for the classical Lorentzian solutions {b{t),x{t)) that are the integral curves of 
^ = -Im(J). 

With this data we then integrate the Lorentzian equations to find the classical solutions 
labeled by {b^,x*)- Later we find the probabilities for these Lorentzian histories. The com- 
plete ensemble of classical predictions is the union of those from all the extremizing solutions 
that contribute to the semiclassical approximation to the no-boundary wave function. 

Choosing a different value of 6=,, to implement this procedure simply means that the same 
Lorentzian solution will be labeled by a different value x* of the x at which it intersects the 
new surface in minisuperspace. 
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FIG. 18: The parameter 7 = lm{(j)Q) /Ile{(pQ) plotted along the vertical contour x = 7r/2,y. This is 
the ratio required to have (p real along this contour. On the left the curves reading from bottom to 
top for equal to .2, .4, .6 and .8. The ratio becomes constant at large y. On the right the single 
value /x = 1.5 is plotted showing the generic lack of stabilization at large y. (cf Figured)) 
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FIG. 19: The first figure shows Lorentzian histories of the scalar field for various values of p,. The 
value of is plotted vertically. Starting from the top on the left and moving downward the 

values of p, are .2, .4, .6 and .8. The solutions are not generally time-symmetric although the 
ensemble of Lorentzian solutions is time-symmetric. The second figure shows the time-asymmetry 
parameter tj (j6.2p as a function of fi. (cf. Figure [T3l) 



Zeroth Order — No scalar Field: The Lorentzian equation for the scale factor b{t) is 

(A8) 



d6 



62 - 1. 



There is only one solution to this equation which is the scale factor for empty de Sitter 
space. Choosing the origin of t to be the time of the bounce this is 



b{t) = cosh(t) 



(A9) 



It is easy to verify that this solution satisfies the zeroth order versions of the Cauchy data 
f l4.16p determined by the complex solution flAip . namely 



dh/dt 



(AlO) 
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where t;* = 7r/2 + iF*, = "K* for one extremizing solution and v^, = it/2 — iY^, = — "K* 
for the other. 

First Order in the Scalar Field: The Lorentzian equation for the scalar field is fl4.13bp . 



1 d (t^dx\ , _2 



where h{t) is the scale factor (]A9p determined in zeroth order. For each contributing ex- 
tremizing solution the Cauchy data along the 6 = 6* surface specified by (14.16bp is 



X(t*) = x*, rfxMli^ = - lm(0) , (A12) 

Let's first consider the extremizing solution defined by v^, = tt/2 + iY^,, t^ = Y^. It 
is straightforward to see that the following is the Lorentizian solution with the boundary 
conditions flA12p . 

X{t) = X* Re[G'(7r/2 + it)/G{n/2 + it,)] (A13) 

provided is identified with Kf. Using flA4bp the Lorentzian solutions can also be 
parametrized by the value of (f)o{oo) as described in the main text for the non-perturbative 
case. Figure [19] shows a few examples of x(^)/0o(oo). 

Another set of extremizing solutions is defined by = 7r/2 — zY"*. This subensemble of 
Lorntzian histories is just the time [t) reversed of the one defined by r;^, = 7r/2 + iY^,. The 
whole ensemble of classical Lorentzian solutions is therefore time-symmetric. The individual 
solututions are generally not as Figure [19] shows. 



3. Perturbative Action 

The action function 7(6, x) determines both when classical Lorentzian histories are pre- 
dicted by the NBWF through fl3.13p . and, if so, what their probabilities are through 
exp{—2Ipi). In perturbation theory the action can be expanded in powers of viz. 

/(6,x) = /(°H6) + /(2)(6,x) + --- (A14) 

where the second term is proportional to x^- The terms I^^\b) and I^'^\h^ x) ^'^^ determined 
by evaluating the action integral (14.80 to quadratic order in x using the perturbative complex 
extremizing solutions found in the first part of this Appendix. The integral in that expression 
is carried out along a contour from r = to an endpoint t = v = X iY corresponding to 
the given value of (6, x)- lu the following we carry out this perturbative evaluation focussing 
exclusively on the range 6 > 1 needed for classical prediction. 

Zeroth Order — No Scalar Field: The extremizing solution is given by flAip . The endpoint 
ist; = 7r/2 + i cosh~^(6). The result of carrying out the integral is 

I^'\h) = -^,[l~^{h'~lf"]. (A15) 

Second Order in the Scalar Field: The action integral (14.81) depends on </'(r), a(r) and 
V. There are perturbations in all three. Only the linear order solution for 0(r) is needed 
to evaluate the field contribution to the action to quadratic order in x- ^ turns out that 
the second order perturbations of the scale factor and of the endpoint cancel essentially as a 
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FIG. 20: Curves of constant real part of the action in minisuperspace. When the conditions for 
classicahty are satisfied these wih be classical Lorentzian solutions for large b. That will be the 
case for the curves on the left with ft < 1, but not those on the right for which /2 > 1. Since the real 
part of the zeroth order action contributes an overall additive constant — 7r/(2i?^) it is convenient 
to label the constant action curves by the value of {IH"^ /37r)I^'^\ [cf lXTT] . Reading from left to 
right the values of (2ff^/37r)/(^) shown are .01, .05, .1, .2, .3, .4, .5 and .6 for the = .5 on the 
left. For p, = 1.5 on the right the values are .01, .1, 1, 10, 100. Note that the 6-axis starts from 1 
which is the value of b at the bounce. 



consequence of reparametrization invariance. Integrating the part of the action by parts 
and using the equation of motion for then gives the following simple expression for I^"^^: 

^^'^ib,x) = ^,b'^iv)^iv) . (A16) 



Explicitly, using (IA4b|l . this is 



I^'Hb y) = i!Ly2^3 g[^/2 + ^cosh-^(&)] ^ i!L^2^3^.^^ 

It is straightforward although laborious to check that this perturbation in the action satisfies 
the Hamilton- Jacobi equation fl3.15ap expanded to second order in x as a consequence of 
the equation of motion fl4.9bp . 

Figure [20] and Figure [21] show curves of constant real part of the action plotted in (6, x) 
and {X, Y) coordinates on minisuperspace for typical values of p, less and greater than 1. 
In regions of superspace where the classicality condition (13.131) is satisfied these curves are 
integral curves of the imaginary part of the action —S{b, x) to a good approximation and 
would represent the predicted classical Lorentzian histories to that approximation. In each 
case the curves display the decay of the scalar field with the expansion of the universe that 
is a property of Lorentzian solutions. 



4. Perturbative Classicality and Perturbative Probabilities 

As discussed in Section IIII B| classical Lorentzian histories are predicted when there is a 
region of minisuperspace where the gradients of the real part of the action x) are all 
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FIG. 21: Curves of constant real part of the action in the complex r-plane. The endpoint values 
{X,Y) provide a set of coordinates for minisuperspace that are alternatives to (b,x)- The unper- 
turbed curve where b is real is the vertical curve at X = 7r/2 along which IH^Ir/Stt = -1/3. 
Perturbations in that curve along which the perturbing matter action is constant are shown here 
for several values of ft. On the left are curves where fi has the values .25, .50, and .75 reading up 
from the lowest to highest. On the right curves corresponding to 1.0 to 1.4 in steps of .1 reading 
left to right. The perturbative calculation of these curves is only valid when they remain close 
to the vertical line at X = 7r/2. For that reason no scale is indicated for the X-axis. For fi < 1 
the curves of constant Ir approach classical solutions. For /x > 1 the oscillation at large y is an 
indication that they do not approach classical solutions. 



small compared with those of minus the imaginary part S{b,x)- A convenient measure of 
this classicality condition fl3.13p is the classicality ratio [cf f l3.13p ] 



C/^(6,x) = |V^/«(6,x)|/|Va5(6,x)|. 



(A18) 



When both these ratios are small the classicality condition fl3.13p is satisfied. In lowest 
non-vanishing perturbation theory order we have 
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(Al9a) 



(A19b) 



The ratio Clf, will be small for small but the ratio Cl^ is independent of x for small x- 

Figure [22] shows these ratios for two values of /i — one below 3/2 and one above. Its 
evident that the classicality condition is not satisfied for Cl^ for the larger value and this is 
true for all values fi > 3/2. By contrast, the condition is satisfied for all /i < 3/2. 

Figure [22] shows the situation with respect to classicality in the yU < 3/2 and fi > 3/2 
regimes in a different way. The real part of the action must become constant along any curve 
in superspace that is a predicted classical history. Were Ir not constant along a Lorentzian 
histories, then exp(— 2/^^) could not be its probability. There is one probability for each 
history. In the approximations that we are using in this paper, the relative probabilities 
of classical Lorentzian histories are given by exp(— 2//?) when the classicality conditions are 
satisfied. For small values of the scalar field this is when /i < 1. The results for the real part 
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FIG. 22: The classicality ratios Cl^ and Cl^ plotted for fi = .75 (top pair) and ^ = 3 (bottom 
pair) in minisuperspace (6, x) where Y = cosh~^(6). In lowest order perturbation theory Clb is 
proportional to x^i the value x = -1 was used for these illustrations. Cl^ is independent of x 
in leading order. The classicality condition is well satisfied for Y ^ 2 when fi = .75. For /i = 3 it 
is not satisfied at all because Cly. never drops much below unity. There are no classical histories 
for /u > 3/2. 



of the action are shown in Fig. [2l]for a range of values of /i < 1. Comparison with Figure 
|5] shows that the perturbation theory is a reasonable approximation for small values of (pQ. 
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FIG. 23: The real part of the action along the classical Lorentzian solution with (pQ = 1 for several 
values of p,. The curves on the left reading bottom top range from // = .3 to = .9 in steps of 
.1. On the right the range from bottom to top is /x = 1.1 to fi = 1.7 in steps of .1. The marked 
qualitative difference between the fi < 1 and ft > 1 is important for classical predictions. On the 
left the real part approaches a constant at large t. Its gradient in this direction is small compared to 
the gradient of S. The curves approach the integral curves for classical, Lorentzian solutions with 
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a constant, its gradient remains comparable to the gradient of S, and, as a consequence, classical 
behavior is not predicted for the scalar field. 
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FIG. 24: The real part of the perturbative Euclidean action Ir (1A14P is plotted for three values 
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